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Abstract — Digital image processing remains a challenging 
domain of programming. All digital images contain some 
degree of noise. Often times this noise is introduced by the 
camera when a picture is taken. Image denoising algorithms 
attempt to remove this noise from the image. In this paper the 
method for image denoising based on the nonlocal means 
(NL-means) algorithm has been implemented and results have 
been developed using matlab coding. The algorithm, called 
nonlocal means (NLM), uses concept of Self-Similarity. Also 
images taken from the digital media like digital camera and the 
image taken from the internet have been compared. The image 
that is taken from the internet has got aligned pixel than the 
image taken from digital media. Experimental results are given 
to demonstrate the superior denoising performance of the 
NL-means denoising technique over various image denoising 
benchmarks. 

Index Terms — ASIC, Image denoising, Non-Local Means 
(NL-means) Algorithm, VHDL. 

I. INTRODUCTION 
Digital Image Processing: - 

Most of the common image processing functions available in 
image analysis systems can be categorized into the following 
four categories: 

1 . Preprocessing 

2. Image Enhancement 

3. Image Transformation 

4. Image Classification and Analysis 

Most denoising algorithms make two assumptions about the 
noisy image. These assumptions can cause blurring and loss 
of detail in the resulting denoised images. The first 
assumption is that the noise contained in the image is white 
noise. This means that the noise contains all frequencies, low 
and high. Because of the higher frequencies, the noise is 
oscillatory or non-smooth. The second assumption is that the 
true image (image without the noise) is smooth or piecewise 
smooth [5]. This means the true image or patches of the true 
image only contain low frequencies. 

Previous methods attempt to separate the image into the 
smooth part (true image) and the oscillatory part (noise) by 
removing the higher frequencies from the lower frequencies. 
However, not all images are smooth. Images can contain fine 
details and structures which have high frequencies. When the 
high frequencies are removed, the high frequency content of 
the true image will be removed along with the high frequency 
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noise because the methods cannot tell the difference between 
the noise and true image [5] [6] .This will result in a loss of 
fine detail in the denoised image. Also, nothing is done to 
remove the low frequency noise from the image. Low 
frequency noise will remain in the image even after denoising. 
Numerous and diverse denoising methods have already been 
proposed in the past decades, just to name a few algorithms: 
total variation[15], bilateral filter or kernel regression[2][16] 
and wavelet-based techniques. [3] [4] [19] [20]. All of these 
methods estimate the denoised pixel value based on the 
information provided in a surrounding local limited window. 
Unlike these local denoising methods, non-local methods 
estimate the noisy pixel is replaced based on the information 
of the whole image. Because of this loss of detail Baudes et al. 
have developed the non-local means algorithm [5] [6] [7]. 

The rest of this paper is organized as follow. In section 2, we 
introduced the non-local means algorithm. Section 3 provides 
experimental work and simulation and section 4 provides 
results and some discussion about above mentioned non-local 
means algorithm. The last section conclude the whole paper. 

Related Work:- 

A. Buades, B. Coll, and J. Morel proposed a non-local 
algorithm for image denoising [7]. The search for efficient 
image denoising methods is still a valid challenge at the 
crossing of functional analysis and statistics. In spite of the 
sophistication of the recently proposed methods, most 
algorithms have not yet attained a desirable level of 
applicability. All show an outstanding performance when the 
image model corresponds to the algorithm assumptions but 
fail in general and create artifacts or remove image fine 
structures. The main focus of this paper is, first, to define a 
general mathematical and experimental methodology to 
compare and classify classical image denoising algorithms 
and, second, to propose a nonlocal means (NL-means) 
algorithm addressing the preservation of structure in a digital 
image. The mathematical analysis is based on the analysis of 
the “method noise,” defined as the difference between a 
digital image and its denoised version. The NL-means 
algorithm is proven to be asymptotically optimal under a 
generic statistical image model. The denoising performance of 
all considered methods are compared in four ways; 
mathematical: asymptotic order of magnitude of the method 
noise under regularity assumptions; perceptual-mathematical: 
the algorithms artifacts and their explanation as a violation of 
the image model; quantitative experimental: by tables of L2 
distances of the denoised version to the original image. The 
most powerful evaluation method seems, however, to be the 
visualization of the method noise on natural images. 

A. Efros and T. Leung[9] proposed Texture synthesis by non 
parametric sampling. Textures can often more easily be 
described as a composition of subtextures than as a single 
texture. The paper proposes a way to model and synthesize 
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such “composite textures”, where the layout of the different 
subtextures is itself modeled as a texture, which can be 
generated automatically. This procedure comprises manual or 
unsupervised texture segmentation to learn the spatial layout 
of the composite texture and the extraction of models for each 
of the subtextures. Synthesis of a composite texture includes 
the generation of a layout texture, which is subsequently filled 
in with the appropriate subtextures. This scheme is refined 
further by also including interactions between neighboring 
subtextures.. 


Image Denoising :- 

Digital images are often contaminated by noise during the 
acquisition. Image denoising aims at attenuating the noise 
while retaining the image content. The topic has been 
intensively studied during the last two decades and numerous 
algorithms have been proposed and lead to brilliant success. 

, N ( 1, if m = n 
(Vm’Vn)- | q, otherwise (1) 

A thresholding estimator project the noisy signal to the basis, 
and reconstructs the de-noised signal with the transform 
coefficients larger than the threshold T : 


/= 

Where 




neighborhoods top’s neighborhood. In the NLM method, the 
denoised value of a pixel is determined by pixels with 
similar neighborhoods. 



Fig 1 : Example of self-similarity in an image. Pixels p and ql 
have similar neighborhoods, but pixels p and q2 do not have 
similar neighborhoods. Because of this, pixel ql will have a 
stronger influence on the denoised value of p than q2 

Non-local Means Method: - 

Each pixel p of the non-local means denoised image is 
computed with the following formula: 

NL{¥){p)= ^w(p,q)V(q) ( 6 ) 


Pt 


M = {ft 


_ JX if |x| > T 
otherwise 


( 3 ) 


is a thresholding operator. 

The mean square error (MSE) of the thresholding estimate can 
be written as 
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The first and seconds terms are respectably the bias and 
variance of the estimate, when the noice is 
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where l{ • }l denotes the cardinal of the set { • } 


II. METHODOLOGY 

In this paper image transformation based on pixel processing 
has been done, which includes image denoising. The 
Self-Similarity concept was originally developed by Efros 
and Leung for texture synthesis [9]. The NLM method 
proposed by Buades [6] [7] is based on the same concept. This 
concept is better explained through an example given in figure 
1. The figure shows three pixels p, ql, and q2 and their 
respective neighborhoods. It can be seen that the 
neighborhoods of pixels p and ql are much more similar than 
the neighborhoods of pixels p and q2 . In fact, to the naked 
eye the neighborhoods of pixels p and q2 do not seem 
to be similar at all. In an image adjacent pixels are most 
likely to have similar neighborhoods. But, if there is a 
structure in the image, non-adjacent pixels will also have 
similar neighborhoods. Figure 1 illustrates this idea clearly. 
Most of the pixels in the same column as p will have similar 


following conditions 0 D w( p, q ) ^ 1 and 



w(p,q) □ 1 . 


Each pixel is a weighted average of all the pixels in the image. 
The weights are based on the similarity between the 
neighborhoods of pixels p and q [1, 2]. For example, in 
Figure 1 above the weight w(p,ql) is much greater than 
w(p,q2) because pixels p and ql have similar neighborhoods 
and pixels p and q2 do not have similar neighborhoods. In 
define let Ni be the square neighborhood centred about pixel 
i with a user-defined radius Rsim. To compute the 
similarity between two neighborhoods take the weighted sum 
of squares difference between the two neighborhoods or as a 
formula 


d {p,q) = || V { N p ) - V ( N* )\\ 2r [1, 2]. (7) 

Here F is the neighbourhood filter applied to the squared 
difference of the neighborhoods and will be further 
discussed later in this section. The weights can then be 
computed using the following formula: 

1 -<Hp y q) 

W(p , q) = e h ' (8) 

Z(/>) 


Z(p) is the normalizing constant defined as 

Z -d(p,q) 

^ h [1 ’ 2] ' 

h is the weight-decay control parameter. 



As previously mentioned, F is the neighborhood filter with 
radius Rsim. The weights of F are computed by the following 
formula: 
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( 10 ) 


where m is the distance the weight is from the center of the 
filter. The filter gives more weight to pixels near the center of 
the neighborhood, and less weight to pixels near the edge of 
the neighborhood. The center weight of F has the same 
weight as the pixels with a distance of one [7]. Despite the 
filter's unique shape, the weights of filter F do sum up to one. 
Equation 1 from above does have a special case when q=p 
This is because the weight w(p, p) can be much greater than 
the weights from every other pixel in the image. By definition 
this makes sense because every neighborhood is similar to 
itself. To prevent pixel p from over- weighing itself let w(p ,p) 
be equal to the maximum weight of the other pixels, or in 
more mathematical terms 


w (p , p) = max w(p , q) \p ± q\ (11) 


III. EXPERIMENTAL WORK AND SIMULATION 

In this section , to verify the characteristics and performance 
of non-local means algorithm, a variety of simulation are 
carried out on the 5 12*5 12 bit gray scale image( place.jpg & 
Einstein.jpg) as shown in fig. 2. All simulations are performed 
in MATLAB 7.0.[17][18] 

In the simulations, firstly images is corrupted by noise and 
then it is denoised using Non-Local Means Algorithm. 

Fig. 2 showed conversion of RGB image into gray scale 
image. 


Original Noisy 



Fig 3 : Image Denoising 


Original Image Noisy Image Noisy image 



Original Image Noisy Image Noisy image 

I" I I I I I '| 
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Fig 4: Original and noisy image in grey scale and noisy 

image in RGB 



Fig 2: Conversion of RGB image to Gray Scale Image 


Finally in fig. 3-6, we present some qualitative results 
obtained by simulation using Matlab 7.0. fig. 3 shows a 
original and noisy image and its denoising. Fig. 4 provides 
original and noisy image in gray scale and noisy image in 
RGB of place and Einstein. 


The aim of fig. 5 is to compare the quality of denoised image 
obtained by various denoising method by visual inspection. 
Fig. 6 shows the blurred image & its restoration using Non- 
Local Means algorithm for different NSR. 


Original image With noise 


Mean filter Avg. 


Wiener 
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Fig 5 : Original, noisy, mean filter, avg , wiener filter image 

of palace and Einstein 

Original Image Simulate Blur and Noise 



Restoration of Blurred. Noisy Restoration of Blurred. Noisy 
Image Using NSR=0 Image Using Estimated NSR 


Original Image Simulate Blur and Noise 


Restoration of Blurred. Noisy Restoration of Blurred. Noisy 
Image Using NSR=0 Image Using Estimated NSR 


I V . RESULTS AND DISCUSSIONS 

There are several transforms that do image denoising. In this 
paper non local means (NL-means) algorithm[5][6][7] is 
used. The denoised images obtained with various 
algorithms are shown in Fig. 5 for visual comparison. Pixel 
based processing is easy to perform as well as it will give 
accurate results in comparison to other methods. Two images 
have been taken in which one is available on system and the 
other which is taken from the digital media and then 
downloaded to the computer. Difference is occurred in the 
processing of two images i.e. the image which is already 
available has got aligned pixels than the image that is 
downloaded from the digital media. 


V. CONCLUSIONS 

This paper gives a generalized method for image denoising. 
Then in depth talk about the non-local means algorithm[7] 
for removing noise from digital image was given. The based 
on simulation results, obtained by Matlab 7.0. Non-local 
means algorithm for image denoising is analyzed on place 
and Einstein images. Also it was shown by experimental 
results that the performance of the NL-means 
algorithm[5][6] outperforms several well-known denoising 
benchmarks in terms of the NSR value. 

In future the matlab code can be converted to VHDL code 
and implemented on FPGA kit in order to develop ASIC 
(application specific IC) for image transformation and 
analysis. ASIC can be made for doing the specific work of 
image denoising ,so a person who don’t know any algorithm 
for image denoising are also capable of doing it. 
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